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Abstract 

We study the freezing kinetics of colloidal polycrystals made of micelles of Pluronic F108, a 
thermosensitive copolymer, to which a small amount of silica nanoparticles of size comparable to 
that of the micelles are added. We use rheology and calorimetry to measure T c , the crystallization 
temperature, and find that T c increases with the heating rate T used to crystallize the sample. To 
rationalize our results, we first use viscosity measurements to establish a linear mapping between 
temperature T and the effective volume fraction, 99, of the micelles, treated as hard spheres. Next, 
we reproduce the experimental T dependence of the crystallization temperature with numerical 
calculations based on standard models for the nucleation and growth of hard spheres crystals, 
classical nucleation theory and the Johnson-Mehl-Avrami-Kolmogorov theory. The models have 
been adapted to account for the peculiarities of our experiments: the presence of nanoparticles 
that are expelled in the grain boundaries, and the steady increase of T and hence (p during the 
experiment. We moreover show that the polycrystal grain size obtained from the calculations 
is in good agreement with light microscopy data. Finally, we find that the (p dependence of the 
nucleation rate for the micellar polycrystal is in remarkable quantitative agreement with that found 
in previous experiments on colloidal hard spheres. These results suggests that deep analogies exist 
between hard-sphere colloidal crystals and Pluronics micellar crystals, in spite of the difference 
in particle softness. More generally, our results demonstrate that crystallization processes can be 
quantitatively probed using standard rheometry. 

PACS numbers: 
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I. INTRODUCTION 



The crystallization of colloidal suspensions has been investigated in detail because of its 
intrinsic interest (e.g. in protein crystallization and as a means to understand crystal- 
lization in atomic and molecular systems, for which colloids are often regarded as a model 
system 2] whose characteristic time and length scales are more readily accessible. Thanks 
to advanced light scattering techniques and scanning confocal microscopy, experiments have 
provided unprecedented information on the nucleation and growth of crystals, including the 
volume fraction dependent nucleation rate or the structure of the nuclei flQ. Most ex- 
periments have focussed on hard-sphere or charged colloids, although more recently more 
complex systems, such as thermosensitive colloidal microgels or mixtures of colloids and 
polymers have emerged 12|-|l6|. In most experiments on colloids, crystallization has been 
investigated following a quench from the fluid phase to the crystal phase. This is usually 



achieved by shear melting the samples and then stopping abruptly the shear 



16|4l8|, or 



by changing abruptly the temperature for mixtures of colloids and thermosensitive surfac- 
tant micelles [0 that induce attractive depletion interactions. To the best of our knowledge, 
a slow transition from the fluid to the crystalline phase has not been investigated yet, in 
spite of the obvious relevance of such a protocol for the comparison with atomic and molec- 
ular materials. In fact, such an experiment would require tuning precisely and in situ the 
particle volume fraction, which is in general a very challenging task. Thermoswellable mi- 
crogel particles could in principle be used 19], but typical swelling ratios are too small 
to allow a wide range of volume fractions to be covered by a sample prepared at a single 



20], where the 



microgel concentration. A different approach has been proposed in Ref. 
slow evaporation of the solvent in a microfluidic chip is exploited. In this paper, we use the 
temperature-dependent formation of spherical micelles in Pluronics block copolymer solu- 
tions to investigate the crystallization kinetics in a colloidal system whose volume fraction 
is increased at a constant rate. 

Triblock copolymers of Pluronics type are a class of commercial polymers made of two 
lateral water-soluble polyethylene oxide blocks and a central polypropylene oxide block, 
whose degree of hydrophobicity can be tuned by varying the temperature, T. Depending 
on polymer concentration and temperature, the Pluronics copolymer can form spherical 
micelles that are sufficiently monodisperse to crystallize at high volume fraction. Due to 
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their amphiphilic character, thermoresponsive properties, and biocompatibility, the Pluron- 
ics copolymer are extensively used in industrial applications as antifoaming or thickener 
agents for instance, and have been recognized as having potentials in biomedical sciences 



21 



22j | . Because these materials exhibit the viscoelastic properties of gels at room tem- 



perature and can flow at lower temperature, they have also been considered as potentially 
interesting hydrogels in electrophoresis experiments 23N26|. Pluronics samples are also at 



the heart of numerous more fundamental investigations. Depending on the lengths of the 
different building blocks, the temperature, and the concentration, Pluronic copolymer can 
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self-assemble and form liquid crystalline and crystalline phases in aqueous solvents 
In particular, crystalline phases made of spherical polymer micelles arranged on cubic lat- 
tices in water are frequently obtained. These phases have been used to template a crystalline 
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30] or as convenient systems to investigate the non-linear rhe- 
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32J. However, their crystallization kinetics 



assembly of nanoparticles 
ology and flow properties of polycrystals 
have been hardly investigated. 

We have recently taken advantage of the thermophysical properties of Pluronics micelles 
to study the effect of the thermal history of Pluronics micellar crystal on the microstructure 



of polycrystalline samples. We have shown using confocal microscopy [33] that the average 
grain size of Pluronics polycrystals can be tuned by changing the heating rate, a property 
well established for molecular and atomic systems, but hardly investigated in colloidal ma- 
terials. In this paper, we deepen the experimental study of the role of the heating rate on 
Pluronics concentrated samples 34|. We use standard rheometry to follow the crystalliza- 
tion processes while T is continuously increased at a fixed rate, showing that varying the 
temperature of the Pluronics sample results in a continuous increase of the volume fraction 
of the micelles. The crystallization temperature is found to increases with T: this behavior 
is rationalized using standard theories for the nucleation and growth of crystals developed 
for atomic and molecular systems and extended to colloidal materials. Our analysis allows 
kinetics crystallization parameters to be accessed and compared to experimental data for 
hard-sphere suspensions. Finally, we provide a quantitative link to the microstructure of 



35]. 



the colloidal polycrystal, a property hardly explored in colloidal crystals 

The paper is organized as follows. We describe the system and the experimental tech- 
niques in Sec. [Ill and present our results in Sec. II III Section [TV] is devoted to a detailed pre- 
sentation of the modeling of the nucleation and growth of crystallites under non-isothermal 
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conditions and of the resulting polycrystal texture. Finally, in Sec. [V] we quantitatively 
analyze our experimental results in the framework of our model and critically compare them 
to numerical and experimental results for colloidal hard-sphere suspensions. 



II. MATERIALS AND METHODS 



A. Samples 

The copolymer polycrystals are composed of an aqueous suspension of Pluronic F108, a 
commercial PEO-PPO-PEO triblock copolymer purchased by Serva Electrophoresis GmbH, 
where PEO and PPO denote polyethylene oxide and polypropylene oxide, respectively. Each 
PEO block is made of 132 monomers, and the central PPO block is made of 52 monomers. 
The copolymer mass fraction is fixed at 34%. The copolymer is fully dissolved at T ~ 0°C; 
upon heating, the PPO central block becomes increasingly hydrophobic, resulting in the 



formation of micelles with a diameter of 22 nm 



36j | . whose number increases with T, even- 



tually leading to a crystalline phase at room temperature due to the packing of the micelles 



37j |. The copolymer polycrystals are seeded with 1 % volume fraction nanoparticles. We 
used Bindzil® plain silica nanoparticles (kind gift from Eka Chemical, sample type 40/130), 
with an average diameter of 30 nm and a relative polydispersity of 19%, as determined by 
transmission electron microcopy. We have recently shown by small angle neutron scatter- 
that the presence of nanoparticle (NPs) does not perturb the crystalline order of 



mg 
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the copolymer micelles. 



B. Experimental techniques 

Rheology measurements have been performed using a stress-controlled rheometer (Phys- 
ica UDS 200) equipped with a Couette cell. A thin layer of low viscosity silicon oil is spread 
on top of the sample to prevent water evaporation. Temperature is controlled by a circulat- 
ing water bath. A temperature ramp from 3 to 23°C at a given rate, T, is imposed to the 
circulating water. The actual sample temperature has been checked for each imposed T by 
inserting a temperature probe directly in the sample confined between the cup and the bob 
of the Couette cell. Typically, T has been varied between 10~ 3 °C min" 1 and 2 °C min -1 . 
For the fastest ramps, T > 0.007 °C min -1 , the same rate is applied in the full range of 
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temperatures, 3 °C < T < 23 °C. In order to reduce the time required to prepare a sam- 
ple, slower ramps are run in three steps: a relatively fast ramp T = 0.5 °C min -1 is first 
imposed, up to T = 12 °C. The desired slow ramp at T < 7 x 10 -3 °C min" 1 is then 
imposed in the intermediate temperature range 12 °C — 17 °C, during which the sample is 
fully solidified. Finally, the temperature is raised to 23° C at a rate T = 7 x 10 -3 °C min -1 . 
We have checked that such a composite ramp provides the same rheological response and 
crystallization temperature as a non-composite ramp with a fixed temperature change rate, 
equal to that imposed in the intermediate temperature range from 12 to 17°C min -1 . 

Complementary calorimetric experiments have been performed on a Micro-DSC III from 
Setaram, using deionized water as a reference sample. Cooling and heating ramps between 
0.5 and 25.7 °C at different imposed speeds, T, have been used. An isothermal step (1 h 
at 0.5°C) is used to thermally equilibrate the heat flow before starting the heating ramp, 
and an isothermal step (900 s at 25.7°C) is used before starting the cooling ramp. In the 
following, we show the data acquired during the heating ramps. 

Sample imaging has been performed with an upright Leica microscope equipped with 
an air x63 objective using differential interference contrast. The samples are confined in 
chambers made by two coverslips separated by a 250 /imm thick 16 x 16 mm 2 double- adhesive 
gene frame (Thermo Scientific). Samples are prepared by placing the sample chamber in a 
copper container immersed in a Haake thermal bath, whose temperature is raised from 3 °C 
to 20 °C at a controlled temperature rate. 
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Figure 1: (Color online) Storage (squares) and loss (circles) moduli as a function of temperature 
during a temperature ramp at a fixed positive rate. In (a) the temperature rate is 0.007 °C min -1 , 
in (b) data are shown for three different temperature rates (0.007, 0.1 and 0.5 °C min -1 ). The 
strain amplitude is 0.1% and the frequency 0.5 Hz. Inset: frequency dependence of the storage 
(squares) and loss (circles) moduli in the fluid phase at low temperature (top) and in the crystalline 
phase at room temperature (bottom). 

III. EXPERIMENTAL RESULTS 

A. Mapping between temperature and volume fraction 

The complex modulus in the linear regime is measured as a function of time at a fixed 
frequency, / = 0.5 Hz, while a temperature ramp is imposed at rate T, from T = 3 °C to 
T = 23 °C (fig. [T^i). At low temperature the sample is fluid-like with a loss modulus, G", 
of the order of 1 Pa and a storage modulus, G', of the order of 0.1 Pa. Above ~ 4 °C, G" 
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increases smoothly while G' is constant. This regime is followed by an abrupt increase of 
both G' and G", and the storage modulus eventually exceeds the loss modulus. At the final 
temperature G' is of the order of 12000 Pa and is at least one order of magnitude larger than 
G" . The very abrupt fluid-to-solid transition is the signature of the sample crystallization. 
We define T c as the temperature at which the storage and the loss modulus cross each 
other, and identify T c as the crystallization temperature. We show in the inset of fig. [T^l 
the frequency dependence of the complex modulus at low (T = 4°C) and room (T = 23°C) 
temperature. At low temperature the sample exhibits a typical terminal behavior of a 
Maxwell fluid, G" = rjoco with an effective viscosity rj = 0.32 Pa.s, angular frequency 
oj = 27r/, and G' = 7] q tu 2 , with a relaxation time r = 0.056 s. At room temperature, 
the behavior is typical of a viscoelastic solid, with a frequency-independent storage modulus 
more than one order of magnitude larger than the weakly frequency-dependent loss modulus. 
Note that the transition is reversible, as the subsequent cooling of the viscoelastic solid leads 
back to the initial low-temperature Maxwell fluid. 

In order to rationalize our experimental results, we propose to use existing theories for 
colloidal suspensions, for which the standard control parameter is the volume fraction of 
colloids. To this aim, the first step is to provide a mapping between temperature and 
volume fraction for the sample investigated here. A mapping has been derived in js^l by 
modeling the structure factor of the micelles with that of hard spheres. For a system very 
similar to ours, the authors found that the volume fraction of micelles varies linearly with 
temperature. Therefore, we assume that the volume fraction of the micelles varies linearly 
with T, 

<f> = a{T-T ), (1) 



and we propose 39[ to determine the parameters To and a by comparing our rheology 
data to viscosity measurements for colloidal hard-sphere suspensions [9|, |40] . The complex 
viscosity, rj = V G' 2 + G" 2 , is calculated and normalized by the complex viscosity at low 
temperature (r/o = 0.32 Pa.s) and plotted as a function of temperature in fig. [2j Up to about 
~ 4.2 °C, the viscosity is constant, and then it increases smoothly up to a temperature 
around 15 °C, where it has grown by one order of magnitude. We find the best collapse of 
the two sets of data for To = 4.2 °C and a = 0.0455 °C _1 (fig. EJ). A reasonable agreement 
is found up to a 5 fold increase of the viscosity corresponding to a volume fraction of about 
0.35. For volume fractions in the range 0.35 — 0.50 (before crystallization), the reduced 
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Figure 2: (Color online) Complex viscosity normalized by the viscosity of the solvent at low 
temperature (0.32 Pa.s) as a function of temperature (bottom axis) and effective volume fraction 
(top axis), for samples prepared with different heating rates as indicated in the legend. Inset: 
Viscosity of the suspension normalized by the viscosity of the solvent versus volume fraction for 



our samples (blue triangles) and for two independent measurements 
suspensions (red symbols). 



40J of colloidal hard spheres 



viscosity of the copolymer sample increases more gently than that of supercooled hard 
spheres suspensions. This discrepancy may have several origins: first of all, as micelles 
form, the "solvent viscosity" might decrease due to polymer consumption. Moreover, the 
aggregation number and the size of the micelles might slightly increase with temperature 41], 
inducing slight changes in the effective volume fraction. Finally, the micelle softness, as 
compared to hard sphere particles, might play a role. In spite of these discrepancies at 
<f) > 0.35, we find that crystallization occurs for an equivalent hard sphere volume fraction 
on the order of 0.52 — 0.56, a range similar to that observed experimentally for hard-sphere 
suspensions. Note that a similar linear mapping between temperature and equivalent hard 
sphere volume fraction was previously derived 37j, by modeling the structure factor of the 
micelles with that of hard spheres. 

The temperature dependence of the complex modulus (fig. [T^,) can be rationalized using 
the mapping thus established: from T = T to T = T c , the sample becomes more and more 
viscous (G" increases) as the volume fraction of micelles increases continuously. The sample 
elasticity is not affected (G' remains small and constant), until crystallization occurs due to 
micelle crowding, yielding the abrupt upward jump of both G' and G" . 
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B. Effect of the heating rate on the crystallization temperature 



Visco-elasticity data measured for different heating rates are shown in figs. [Tb and [2j 
Prior to crystallization, all data collapse nicely, indicating that the rate of increase of the 
micelle volume fraction is not controlled by kinetics factors but is instead controlled by 
thermodynamics and thus depends only on temperature. Similarly, data collapse in the 
solid crystalline phase at high temperature/ volume fraction, showing that the mechanical 
properties of the solid phase are robust with respect to a change of the solidification protocol. 
Note however that there are slight differences in the loss modulus, which will be analyzed 
in details elsewhere. By contrast, we find that the crystallization temperature depends 
significantly on T. Although this is not a priori surprising since crystallization involves 
kinetic processes, it is interesting to observe that simple rheology experiments can capture 
these effects. Figures [Tb and [2] show that the crystallization temperature is shifted toward 
higher temperatures as T increases. The evolution with the heating rate of the crystallization 
temperature, T c , as measured by rheology, is plotted in fig. [3j for T spanning more than 
three orders of magnitude. Numerical values are in the range (15.5 — 19.5) °C. Although 
rather scattered, data clearly indicate a decrease of T c when T decreases leading eventually 
to an "equilibrium" value T c ~ 15.5 °C for very slow heating rates, essentially independent 
of t. 

We complement the rheology data with DSC measurements of the crystallization tem- 
perature. A complete DSC curve is shown in fig. HI One can distinguish a very broad and 
deep inverse peak, due to micellization, which is centered around T = 5 °C. This peak is 



followed 
micelles 



39 



a low amplitude peak around 16 °C, which is due to the crystallization of the 



42]. Differential scanning calorimetry data acquired at different T, although 



in a range smaller than for rheology experiments due to experimental limitations, show a 
similar shift toward higher temperature of crystallization as T increases. Numerical values 
of the crystallization temperature extracted from DSC are plotted together with the rheol- 
ogy data in fig. El Data measured with the two techniques are compatible, supporting the 
method used to extract a crystallization temperature from visco-elasticity measurements. 
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Figure 3: (Color online) Crystallization temperature as a function of the heating rate as measured 
by linear rheology (blue circles) and DSC (green stars). The red line is a fit to the experimental 
data using the model described in the text. 

IV. MODELING THE NUCLEATION AND GROWTH OF CRYSTALS 



The fact that crystallization occurs at a higher temperature (higher 0) when solidification 
proceeds faster can be understood by recalling that crystallization is a kinetic process. As 
shown in a so-called time-temperature-transformation diagram 43| (TTT diagram), crys- 
tallization occurs at small supercooling (smaller volume fraction, and hence lower T in our 
case) for a very slow temperature ramp, while it occurs at higher supercooling (smaller vol- 
ume fraction, higher T) for very fast ramps. To quantitatively account for our experimental 
observations, however, one needs to take into account not only the kinetically controlled 
nucleation process but also the growth process of the crystalline phase. We use here stan- 
dard theories for the nucleation and growth of crystalline phases, initially developed for 
atomic and molecular systems and subsequently applied to colloidal suspensions. Moreover, 
we adapt these theories in order to take into account the presence of nanoparticles and the 
non-isothermal conditions under which the sample solidifies. 



A. Crystallization temperature and average crystallite size 



We first recall how an average grain size can be computed, using the Johnson-Mehl- 
Avrami-Kolmogorov (JMAK) theory 4J-|48[. In this theory one assumes that nucleation 
occurs randomly and homogeneously and that the growth speed of the crystallites, Vq, does 
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Figure 4: (Color online) Heat flow versus temperature when temperature is varied from 0.5 to 
25 °C at a rate of 0.1 °C/min. Inset: zoom of the data in the zone when the crystallization peak 
occurs, for data acquired at different temperature rates as indicated in the legend. Data in the 
inset have been arbitrarily shifted along the y axis. ** and * point respectively to the peaks due 
to micellization and crystallization. 



not depend on the extent of the crystallization process and is isotropic. Following JMAK, 
the extended volume fraction of the sample, X e (t), that has crystallized at time t, reads: 



X e {t) 



Air 



Kt) 



n 3 



v G (r')dr' 



dr 



ill 

3 



/> G (r')dr' 



(2) 



is the 



where J(r) is the nucleation rate per unit volume at time r, and 
volume at time t of a grain that has nucleated at time r and has grown from time r to time 
t, with a time-dependent growth velocity vg{t')- 

The extended volume fraction does not take into account the impingement between crys- 
tallites. As shown by Kolmogorov and Avrami, the actual volume fraction of the sample 
that has crystallized is related to X e (t) by 



X(t) 



exp[-X e (t)] 



(3) 



The average grain size R at the end of the crystallization process can be derived from 



the final grain density 



49 
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(4) 
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Figure 5: (Color online) Evolution of the volume fraction of sample that has crystallized, X, with 
temperature (bottom x-axis) and with volume fraction (top x-axis). X is computed with the same 
parameters as the ones used to fit the experimental data (fig. E}. Data are shown for several heating 
rates, as indicated in the caption. 

Here I a is the actual nucleation rate which takes into account the fact that nucleation can 
only occur in regions which have not crystallized yet: 



Ut) 



[1-X(t)]/(t). 



(5) 



The above equations are very general and hold for time-dependent nucleation and growth 
rates, as is the case in our experiments where a continuous heating is imposed. Note that, 
using these equations, the standard result obtained in the case of time-independent nucle- 



ation and growth rates can be recovered [35|. In this case, X e (t) = -y-Ty^t 4 . If one moreover 



neglects impingement and considers the extended volume fraction of crystallized material as 
the actual volume fraction, the time t* at which crystallization stops reads t* ~ (Jt^) -1 / 4 . 
Hence the average crystallite size (eq. HJ) reads R ~ I ^ J* drj ~ (^f) 1 ^, as previously 



35). 



derived but not experimentally verified for colloidal samples 

For a time- varying process, one can use Eqs. (|2J) and to calculate the temporal evolu- 
tion of the crystallinity, provided that expressions for vq and I are given. Before discussing 
in detail how to model the T- (or, equivalently, 0-) dependence of the nucleation and growth 
rates, we anticipate the results for the calculated X(t) in the case of our samples. The 
evolution of the actual volume fraction of the sample that has crystallized is shown as a 
function of temperature (bottom axis) or volume fraction (top axis) in fig. [51 for several T 
spanning three orders of magnitude. We take as time t — the time at which the sample 
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temperature reaches the equilibrium solidification temperature (T = 15.5 °C, <fi = 0.514). 
All curves have sigmoidal shapes and reach asymptotically 1. They are shifted towards 
higher temperatures/volume fractions when T increases, as observed experimentally. Al- 
though X(t) reaches 1 only asymptotically, the curves shown in fig. [5] can be used to define 
a crystallization time, t*, by setting a threshold X t ^ such that X(t*) = X^. Knowing the 
temperature history imposed to the sample, the evaluation of temperature T* at which crys- 
tallization is completed is then straightforward. Having established how the crystallization 
temperature can be in principle calculated, we now need to provide analytical expressions 
for the nucleation and growth rates. 



B. Nucleation rate, growth rate, and the effect of nanoparticles 

We assume that crystallites grow into spherical particles at a rate vq, which is given by 



the Wilson-Frenkel equation 



35|: 



« = ^[l-exp(-|AAi|/* B T)]. (6) 



Here D s is the self diffusion coefficient of the particles (atoms, molecules or colloids) in 
the fluid phase, A is the typical distance over which a particle diffuses to become part of 
a growing crystallite, A/x is the difference in chemical potential between the solid and the 
fluid phases, and ks is Boltzmann's constant. Note that A is expected to be of the order of 
the particle size. All quantities depend on temperature for atomic and molecular systems, 
and on volume fraction for colloidal systems. 

We use the classical nucleation theory (CNT) to evaluate the nucleation rate. Classical 
nucleation theory predicts that the nucleation rate per unit volume, J, is an Arrhenius-like 



function 



5l|: 



I = Texp[-AG*/k B T], (7) 
where AG* is the nucleation barrier and T the kinetic factor 



35|: 
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In Eq. (jHJ), n c is the number of particles in a critical nucleus and pi is the particle num- 
ber density in the fluid phase. In CNT, the critical radius and the nucleation barrier are 
calculated by balancing the surface Gibbs free energy of a crystalline nucleus of radius r, 
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E$ = 47rr 2 7 and the volume Gibbs free energy term, Ey = — -^ps\ A/i|r 3 , where 70 is the 
surface energy between the fluid and the crystal phase and p$ is the number density of the 
particles in the solid phase. Since E a grows as r 2 while Ey decrease as — r 3 , the net change 
of free energy for a nucleus as compared to the fluid phase goes through a maximum: the 
free energy barrier to be overcome for a nucleus to be stable is then AG* = ^ ( ps |A/j|) 2 • 

We now provide a model to account for the presence of a small amount of NPs of diameter 
<Tp that are added to the copolymer suspension, at a volume fraction 0p. We have previously 
observed that upon solidification the NPs concentrate in the grain-boundaries and influence 
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52|. We assume 



the microstructure of the polycrystals 33J. In the following, we neglect the effect of the 
diffusion of NPs on crystal growth, which could occur because of NPs partitioning, and 
thus assume that vq remains unchanged (Eq. By contrast, we argue that the NPs 

influence the barrier for nucleation, AG*, because they are trapped at the interface between 
the crystalline phase and the fluid phase, as observed experimentally 
that the net effect of this trapping is a decrease of the surface free energy of the crystalline 
nuclei. To model this phenomenon, we adapt the CNT calculation of the free energy barrier 
by modifying the surface free energy term Eg- We propose that the surface free energy is 
reduced by an amount E' s proportional to the surface occupied by the NPs trapped at the 
interface. Since these particles were initially in the volume occupied by the crystallite (the 
number of particles in a crystallite of radius r is 8r 3 <j)p/ap), their number scales as r 3 <j)p/ap. 
Thus, E' s ~ apr 3 <j)p/ap, assuming that each NP at the interface contributes to a reduction 
of surface energy proportional to the area it occupied at the interface, fop, so that the 
modified surface free energy reads: 

E s = 4vr 7 or 2 -E' s = 4n lo r 2 - 7o ^^r 3 (9) 

CTp 

where Ep is a proportionality constant accounting for both the partitioning of the NPs (i.e. 
for what fraction of them is actually rejected at the interface of the growing crystallite) and 
the reduction of surface tension due to one single NP. Note that the fact that impurities 
could affect the liquid-solid surface tension was already mentioned in Ref. js3| in the context 



of metallic alloys and in Ref. [54] in the context of phospholipid monolayers. 

By balancing the surface and volume terms, we find a modified equation for the nucleation 
barrier: 

AG * = 16tt To riQ x 

3 HA/i| + §7o£p0pA7 P ] 2 • 
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Interestingly, the NP contribution, although physically due to a change of the surface term, 
can formally be incorporated in the volume term Ey, because it scales with the volume of 
the nucleus. The net effect of the NPs is thus to lower the nucleation barrier thanks to an 
effective increase of the volume term. 



V. ANALYSIS OF THE EXPERIMENTAL DATA 
A. Numerical values of the parameters 

In this subsection, we detail the numerical values and the fitting parameters used in 
evaluating vq and / as defined in Sec. IIVBI Several parameters can be evaluated from 
previous numerical and experimental works on hard-sphere colloidal suspensions, finally 
leaving four fitting parameters. In all calculations, we use the mapping between volume 
fraction and temperature determined experimentally, Eq. (CQ), assuming that this mapping 
is valid in the whole range of volume fractions considered. The number density in the liquid 
phase is pi = with a = 22 nm the micelle diameter. To reduce the number of fitting 
parameters with no significant loss of generality, we assume a constant volume fraction in 
the solid phase (0 S = 0.5), so that ps = ^f- We furthermore fix n c = 100, since for hard 
spheres the number of particles in a critical nucleus has been found numerically to depend 



weakly on volume fraction [35]. Note that the exact value of n c does not have a strong 
influence on the calculations, since only the kinetic parameter T depends on n c , but very 
weakly, T ~ (Eq. [7]). Finally, the micelle diffusion coefficient D s is strongly dependent 
on volume fraction, and hence on temperature. D s could be in principle estimated from the 
viscosity data of Fig. |2j However, viscosity is measured macroscopically on samples where 
crystalline and fluid phases may coexist, while Eqs. (EJ) and ([8]) require the microscopic 
diffusion coefficient in the fluid phase alone. We thus estimate D s from measurements 
of the microscopic structural relaxation time, r a , of supercooled suspensions of colloidal 
hard spheres, assuming D s ~ r" 1 . By fitting the data of Ref. 



law D s = D 



[1 - 2.50 + 1.360 2 )(1 - 1) + 144.570 2 (1 - £ 



2.5 



we find the empirical 
with C = 0.598 and D = 



6 x 10~ 14 m 2 s _1 evaluated from the Stokes-Einstein relation, D s = ^ ^ , using rj = 0.32Pa s. 
This empirical law holds to a very good approximation up to <fi ~ 0.55 and is consistent with 
our viscosity data for the micellar suspension up to the onset of crystallization. Note that 



1(3 



D s vanishes for C = 0.598, which implies, according to Eq. (pQ), that the micelle dynamics 
is totally frozen for T > 17.34 °C. 

The remaining parameters appearing in vq and / are treated as fitting parameters. Ac- 



cording to numerical work on hard spheres 



351 ]. the volume fraction dependence of the 



chemical potential difference between the solid and the liquid can be well adjusted by the 
simple functional form A/x = A A1 /csT(0 cr ^ — 0), with 4> cr ys = 0.49776 and = 14.7. Here 
we consider the prefactor as a fitting parameter and take (f) cr ys = 0.514. With this choice, 
A/x = for T = 15.5 °C, which is the experimental numerical value of T c measured for ex- 
tremely slow heating rates (see fig. [3]), for which the crystallization temperature is nearly 
independent of the heating rate. The surface tension between the fluid and solid phases 
in colloidal suspension has been measured experimentally by several groups. They all find 

In the following, we use the same expression tak- 
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7o = ksT^- with A 7 of order one 
ing A 1 as a fitting parameter. Numerical simulations and experiments show that the typical 
distance over which diffusion occurs is proportional to the colloid size: A = A\a, although 
there is a significant discrepancy between simulation and experiments: A\ is found in the 
range 0.3 — 0.5 in simulations 35 1. while it ranges from 3 to 17 in experiments |6|. In the 
following, we assume A = A\<y with A\ a fitting parameter. The fourth adjusting parameter 
is E p , the proportionality constant that accounts for the effect of NPs in the reduction of 
the surface tension already discussed in Sec. HVBt see Eq. (jTUJ) . 



B. Fit of the experimental data 

We use the model developed above to fit the experimental crystallization temperature 
as a function of temperature increase rate (fig. [3]). We remind that rheology measurements 
yield T c , defined as the temperature at which the sample response becomes predominantly 
solid-like (the storage modulus G' becomes larger than the loss modulus G"). In the model, 
one needs to define a threshold for the crystallinity X in order to define the crystallization 
temperature, T*, since the sample is fully crystallized only for t — > oo (fig. [5]). We identify T c 
with T* choosing a threshold X t ^ = 0.9 and check that the results are essentially threshold- 
independent for 0.5 < Xth < 0.98. The best fit to the experiments is shown as a continuous 
line in fig. [31 showing a good agreement with the data. This adjustment was obtained with 
the following values of the fitting parameters: A^ = 10 ± 2, A 1 = 1.5 ± 0.1, A\ = 0.11 ± 0.3, 



17 



Figure 6: Differential Interference Contrast light microscopy of the samples prepared with different 
heating rates as indicated in the legend. Scale bars: 20/mi 



Ep = 113 ± 5. The values for A^, A^ and Ax are in agreement with those expected from 
numerical simulations and experiments on hard-sphere colloidal suspensions (see Sec. M) . 
They are also reasonably close to the values extracted using the same model for reproducing 
the crystallite size as a function of NP content and temperature rate {33I [A^ = 17 ± 3, 
A 7 = 0.75±0.2, A\ = 0.14±0.05, E P = 126±24], albeit with a different kind of nanoparticle 
impurities. 

As a further check of the soundness of our approach, we compare the T dependence of 
the average crystallite size obtained from the model to that measured experimentally by 
analyzing images of the polycrystalline microstructure (fig. EJ). The images are taken by 
differential interference contrast optical microscopy and the microstructure is visible thanks 
to the presence of nanoparticles that segregate in the grain-boundaries, thereby providing 
the required optical contrast between adjacent grains. To calculate the average grain size, 
R, as a function of T, we use Eq. (j3J) with the same set of fitting parameters as determined 
by fitting T C (T). Experimental and numerical results are shown in fig. [7] as symbols and a 
line, respectively. For T in the range 10~ 4 — 1 °C min , the model predicts R to be in the 
micrometer range and to continuously decrease as T increases, as observed previously 33]. 
Note that the drop of R is significant at relatively fast heating rates, while the grain size 



levels off at about 10 /im for very slow rates. Indeed, as discussed in Ref. [33j, for slow ramps 
the crystallite size is only limited by the presence of nanoparticles. The model is in very 
good agreement with the experimental values of R for samples prepared with temperature 
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0.01 0.1 

T ("C/min) 



Figure 7: (Color online) Average size of the crystallites as a function of the heating rate. The 
symbols correspond to the experimental data and the line is the numerical values calculated with 
the model detailed in the text. The adjusting parameters are the ones extracted from the fit of the 
experimental data of fig. [3j 

rates in the range 0.001 — 0.2 °C min -1 , capturing both the decreasing trend of R and its 
absolute value, thus providing additional support to our analysis. Interestingly, the grain 
size can be measured easily for the lowest heating rates, for which the network of grain 
boundaries is clearly visible in the microscope images (fig. [6^, b). For T = 0.2 °C min -1 
(fig. |6b) the contrast is poorer, but the grain boundaries are still visible. By contrast, the 
sample prepared with the fastest heating rate, T = 2 °C min" 1 (fig. |6ji), appears uniform, 
presumably due to the very small grain size, as predicted by the model, and to a less 
efficient partitioning; of the NP between the bulk and the grain boundaries, as also observed 



in Ref. 



33). 



C. Discussion 



We have used a thermosensitive micellar system, for which the volume micelle volume frac- 
tion can be tuned with temperature, to investigate the crystallization dynamics of colloidal 
suspensions under time-varying volume fraction conditions. The crystallization temperature 
T c was measured using standard rheometry and its dependence on the rate of the tempera- 
ture ramp used to solidify the sample was investigated. The experimental data can be very 
well accounted for by using a standard model for the nucleation and growth of colloidal 
crystals that has been adapted to the case of time- varying volume fraction conditions and to 
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the presence of NP impurities. Our model predicts also how the microstructure of the poly- 
crystal evolves with T, which is found in good quantitative agreement with the experiments, 
thus providing a successful cross-check of our approach. We have therefore demonstrated 
that rheometry experiments combined with standard models for the nucleation and growth 
of colloidal crystals allow quantitative parameters related to the crystallization processes to 
be derived. 

Our analysis was performed assuming that the behavior of the suspension of copolymer 
micelles can be mapped to that of hard-sphere suspensions. To proceed a step further along 
this analogy, we compute the nucleation rate, as given by Eqs. and ffTU]) . using the set 
of adjusting parameters determined in Sec. |Vj and compare it to results from experiments 
and simulations of hard spheres. Data taken from the literature are collected and plotted as 
a function of volume fraction in fig. [HJ together with our model. In order to compare data 
obtained with different particles and solvent viscosity, we plot the normalized nucleation rate 
I* = Jo~ [Dp where a is the colloid diameter and Do its diffusion coefficient in the dilute 
regime 35J. Interestingly, one finds a good agreement between the experimental data of 
hard sphere suspensions and the curve extracted from our model, the agreement being much 
less good with results from numerical simulations of hard sphere suspensions. More quan- 
titatively, our data are slightly shifted towards higher volume fraction as compared to data 
for hard-spheres. We indeed find that the maximum of I* occurs at <ft = 0.58 for copolymer 
micelles, to be compared to <ft — 0.56 for hard-sphere colloids. By contrast, the height of the 
maximum is of the same order of magnitude: we calculate a maximum normalized nucleation 
rate of ~ 2 10~ 8 for our micelles, whereas experiments on hard-spheres yield values in the 
range 2 10~ 5 — 710 -8 . Despite small discrepancies, the remarkable agreement found here 
hints at profound analogies between hard-sphere suspensions and suspensions of Pluronics 
micelles, in spite of the different nature of the materials. In particular, the interaction po- 
tential between micelles is expected to be much softer than that between hard-spheres. The 
softness of the potential is known to have a profound influence on the suspension properties 
near jamming. However, this appears not to be the case when the packing fraction is still 
well below random close packing, as at the onset of crystallization, presumably because the 
particles hardly experience any direct contact. 

To conclude, thanks to our approach combining experiments and modeling, we have 
pointed out analogies between hard-sphere colloidal crystals and Pluronics micellar crystals, 
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Figure 8: (Color online) Normalized nucleation rate, /* = la 5 /Do, as a function of particle volume 

Q Q 

fraction. Symbols are experimental data for hard spheres suspensions of Ref. [7|] (squares), Ref. [3( 
(circles), Ref. [q| (diamonds), Ref. 4] (triangles). Line and symbols are numerical simulations of 



monodisperse (circles) and 5% polydisperse (squares) 



The thick line is our model. Figure adapted from 



rard sphere suspensions, taken from Ref 



in spite of the difference in particle softness. This observation should be of particular interest 
to the numerous groups currently working on soft, thermoswallable microgel particles. In 
addition, we have provided a general framework to rationalize the effect of temperature 
history or, more generally, volume fraction history, on colloidal crystallization. These results 
should be relevant for colloidal systems where volume fraction can be tuned, and in particular 
for the vastly investigated Pluronics gels. Finally, we have demonstrated that crystallization 
processes can be quantitatively probed using standard rheometry, an experimental method 
more readily accessible than most microscopic characterization tools, such as neutron or 
X-ray scattering. 
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